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A discrete version of the two-dimensional inverse scattering prob- 
| lem is considered. On this basis, algebraic transformations for the two- 

' dimensional finite-difference Schrodinger equation are elaborated. General- 

ization of the technique of one-dimensional Darboux transformations for the 
two-dimensional finite-difference Schrodinger equation is presented. Ideas 
of Bargmann-Darboux transformations for the differential one-dimensional 
multichannel Schrodinger equation are used for the two-dimensional lat- 
tice Schrodinger equation. Analytic relationships are established between 
different discrete potentials and the corresponding solutions. 
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I. INTRODUCTION 



> 

in 
O . 

Study of multidimensional and multi-particle objects is qualitatively more complicated 
Q\ • than that of the one-dimensional case. When there is no symmetry in the potentials of 
interaction between particles, these systems are described by partial differential equa- 
tions with unseparable variables. Just for this reason the inverse scattering problem in 
2- and 3-dimensional spaces was formulated by Faddeev and Newton 0|§, Novikov 
and Henkin [|J at a much later than the 1-dimensional problem. The foundation for de- 
veloping the 2-dimensional inverse problem in finite differences was laid by Berezanskii 
H who developed the theory of orthogonal polynomials for the Jacobi infinite matrix. 
The Darboux transformations in quantum mechanics for the Schrodinger differential and 
difference 1-dimensional equation have much in common with the spectral transforma- 



tions for orthogonal polynomials ||. It is therefore expedient to analyze the spectral 
inverse problem for the discrete 2-dimensional Schrodinger equation on the basis of the 
technique of orthogonalization of polynomials. The formulae of the inverse problem with 
degenerate kernels are closely related to Bargmann-Darboux transformations J7|- 0. The 
one-dimensional Darboux transformations has already found wide applications in quan- 
tum mechanics and in the theory of nonlinear integrable systems pTO|JTT] j^JT^] . 

In the present paper, a discrete version of the Gelfand-Levitan inverse spectral problem 
is considered for the two-dimensional lattice Schrodinger equation. The two-dimensional 
finite-difference inverse problem, based on the procedure of orthogonalization of poly- 
nomial vectors, is a generalization of the one-dimensional procedure given in |TJ|. On 



the basis of the obtained formulae of the inverse problem, the relations are derived for 
discrete Bargmann-Darboux transformations in two dimensions. The suggested algebraic 
approach allows one to construct families of discrete potentials in an explicit form and 
the corresponding solutions. 
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II. INVERSE PROBLEM 



Consider the Schrodinger equation whose Hamiltonian is tridiagonal in a certain basis 
with respect to both the coordinate variables n and m 



(Hip) nm = a nm tp(n - 1, m) + a n+lm ij(n + 1, m) 

+b nm ^(n, m - 1) + b nm+ ii/;(n, m + 1) + c nm ${n, m) = \i/;(n, m). 



The coefficients a nm ,6 nm ,c nm 



are assumed to be real and represent discrete potentials; 
•0(n, m) are discrete wave functions, (n,m) is an integer point of the half-plane, n = 
0, 1, 2, .., m = —1, 0, 1, A is a spectral parameter. The index n can vary from to 
oo and the index m can vary from — oo to oo or from to oo; when < n < N and 
< m < M, it is a special class of restricted problems. The equation in finite differences 
(pi) can be represented by the expression with operator coefficients 



= A n m{n - 1) + Vr$>{n) + A n+1 m{n + 1) = A* (n), 



(2) 



when we treat one of the variables, n, as the only discrete coordinate; and the other, m, 
as the channel index. If we set \&(— 1) = 0, the action of the Schrodinger discrete operator 



H on the vector ^ = {#(0), *(1), *(2), ...*(n) 
Jacobi block matrix J on \1/ 



.} is represented by the action of the 
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whose elements V„ and A n are matrices at each fixed n, and each element \I/(n) of the 
vector \l/ corresponds to the vector {ip m (n)}, ipm{n) = m )> i n the other space variable 



"m" 
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It is clear that the operator if is Hermitian. Note that the coupling of equations (|l|) 
or (Q) with respect to m occurs only between immediate neighbours and for this reason 
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the symmetric coupling matrix Ki;(mm') = Vmm'{ n ) m © is tridiagonal and relates the 
functions ip{n,m) at m,m ± 1, unlike the conventional multichannel case. 

From the Jacobi block matrix (0) tridiagonal in the variable n it is seen that every 
vector of the solutions \I/(n) is connected with the vectors of solutions \I/(n ± 1) at the 
neighboring points n ± 1. As a result, if we take nonhomogeneous boundary conditions 
at one end of the interval < n < N, we can obtain the solutions on the whole interval 
by moving by subsequent steps from that end. It turns out that these vectors of solutions 
as functions of the spectral parameter A are polynomials of A with matrix coefficients || 
which can be orthogonalized in the spectral measure. 

The spectral inverse problem is reduced to the construction of potential matrices V n , 
A n and an unknown system of orthonormal polynomials with the use of the known system 
of orthonormal polynomials corresponding to the finite-difference equation (0) but with 

o o 

the known matrices Vn and A n - 



A. Orthogonalization of polynomials 

o 

Introduce auxiliary solutions <p ms (n) = (p s (\, n, m) and (f ms (n) to eq.(hj) with the help 
of the boundary conditions 

</w(-l) = 0; (p ms (0) = S ms ; (5) 

^ms ( 1) 0, ^ms (0) $msi 

o 

where s is a point on the m axis. The solutions ^P s satisfy the same equation ([!]) with the 

o o o o o o o 

known potentials a nm ,c nm =Vmm (n) , b nm +i=V mm+i (n), b nm =Vmm-i {n) 

m+1 

o o o ° x — > ° ° 

a nm Vms (A, n - 1)+ a n+lm ¥ ms (A, n + 1) + Vmm" (n) ¥ m » s (A, n) 

m"=m— 1 

= \¥ ms (\,n). (6) 

As a first boundary condition, we can change the zeroth condition to a more general 
one corresponding to the homogeneous boundary conditions (ip ms (—l) — y? ms (0))/A = 
D s (p ms (0). In what follows we assume the step of finite-difference differentiation A to 

o 

equal 1. According to the conditions the functions <p ms (n) and f ms (n) are equal 
to zero on the lines (— l,m), (0,m) except for the points with m = s on the line n = 

o o 

where y2 ss (n) =<P SS (n) = 1. From these values of (p s and <£ s one can find the functions 

o 

(Pms(n) and ^ ms (n) at all other points of the coordinate network by using the recurrence 
equation ([!]). Since the matrix equation (|2|) is tridiagonal in the variable n, and owing 
to the boundary conditions (^), the vectors $ s (n) = {...<fo s (n),(fi s (n), ...<f ms (n)...} of 
solutions $(n) = ($o(^), ...$ s (n), ..) in the spectral variable A are polynomials of 

the nth degree. At n = 1 the vector $ s (l) is a polynomial of the first degree; at n = 2, the 
vector $ s (2) is a polynomial in A of the second degree, and so on. It appears, however, 
that different elements (p ms (n, A) of the vector of solutions $ s (n) have different degrees 
of polynomials. Because of the three-point coupling with respect to m the functions 
ip ms (n, A) are polynomials in A of the degree n — \s — m\. The maximum nth degree 
belongs to the elements with m = s and the elements <f) ms (n, A) vanish at all s lying out 
of the values (n — m, n + m). The distribution of polynomial degrees is shown in Fig.l. 
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FIG. 1. Nonzeroth 
</> ms (n,A) are inside and on 
boundaries of the cone with 
its vertex at (0,s); the num- 
bers at the nodes of the net 
indicate the degrees of ip ms (n) 
polynomials. 
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It is clear that the vectors of solutions, 
whose elements are polynomials in A 
can be orthogonalized. Berezanskii has 
shown that if the elements a nm , b nm 
and c nm in (|l|) are real and a nm > 
{n = 1, 2, m — ... — 1,0,1,...), one can 
introduce polynomials of the first kind 
(ip ms (X,n) = P s; („, m )(A)) which are or- 
thogonal with the weight of the spectral 
matrix p(A) = {p ss >(\)} 



II 



(7) 



u 

and similarly for V 9 ™,* (A,n) 

_ Co o o 

2^ / Vms {\n)dp ss i (A) <£Vm' (A,n ) = <W<W'- (8) 

ss' 

o 

The spectral matrix elements p ss i (or P ss i) are determined by the boundary values of 

o 

the special solutions i/j(n,m) (or ip (n,m)) obeying the zeroth condition ip(— l,m) = 
i/; m (—l) = (m = .. — 1,0,1...). For example, for p states of the discrete spectrum, 
the matrix p(X) is the sum of p terms formed from productions of the column vectors 
r(A l/ ) = {^s{K)} and the row vectors T^(X V ) = ( r r s (X u )) 

p 

Pss'(A) = ^20(\- A v )7 s (A y )7 S '(A„), 

v=l 

where the elements 7 s (A iy ) = ip s (0, X u ) are defined by ip{n = 0, m) and 9(X — A^) is the 
Heaviside step function equal to 1 at A = \ v and to zero when A ^ \ v . A solution ip(n, m) 
to eq.(P with the zeroth condition (or a more general homogeneous condition )[] can be 
obtained by multiplying the matrix $(A,n) by the vector T(A) = ^(0, A): 



V>(n,m,A)= Prnsin, \)ip s (0, A). 



(9) 



At every fixed n and m the summation over s is finite owing to (p ms (n, A) = outside the 
interval (n — m,n + m). 



* Two solutions of the same equation of second order differ from each other by a normalization 
factor at those points A where both of them exist when one of the boundary conditions is the 
same. 
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B. A discrete version of the Gelfand Levitan inverse 

problem 



Using the procedure of orthogonalization of polynomials we construct unknown poly- 
nomial solutions (p ms (\,n) normalized with the spectral weight p ss i(X) as a linear com- 

o 

bination of the known polynomial solutions ¥> ms (X,n), orthogonal with respect to the 

o 

measure P ss i (A) 

n m+(n—n') 

ip ms (\,n) = H K(n,m;n',m') V m 's (X,n'). (10) 

n'=0 m'=m—(n—n') 

As one can see above, the validity of this relation is a consequence of the fact that both 
the functions, <*P ms (X,n) and ip ms (X,n), obey the same discrete Schrodinger equation ([I]) 
and the same boundary conditions (|5|), given on the lines n — — 1 and n = 0. Therefore, 

o 

the solutions y? ms (A, n) and ^ ms (A, n) are polynomials in A of the same degree n — \s — m\, 
but with different coefficients. The matrix of the polynomial solutions $(A,n), which is 
of degree n, is orthogonal to every polynomial matrix of degree lower than n and hence 

o 

to every $ (A, n 1 ) for n' <n 

J Vms(X,n)dp ss/ (X) Vs'm' (AX) = 0. (11) 

ss' 

Equation (pTDj) is a discrete analog of the Volterra integral equation, in which the coef- 
ficients K(n,m;n',m') are determined by the condition of orthogonality of the vector- 
functions $(A,n) = (...$i(A, n), $ 2 (A, n), $ S (A, n), ...) orthogonal with the spectral 

o o 

measure p(A) to the functions $ (X,n') orthogonal with the weight matrix p (A), when 
n' < n 



Yl J Prns(X,n)(dp ssl (A) - dp ss >{X)) £ s / m / {X,n') = K(n,m;n',m'). 



12) 



The Volterra equations fllCf ) have a triangular form, K{n } m; n\ m') = for n' > n. It is 

o 

easy to see from ( |TD| ) and orthogonality of ^P sm (X,n) (g) that for n' < n 
K(n,m;n',m') =H/ l fms(X,n)d °P SS > (A) <£V m / (X,n'). 

ss' 

Inserting (|T0D into (12) for n' < n we obtain the following system of equations for the 
orthogonalization coefficients K(n, m; n', m') 

K(n, m; n' , m') + K(n, m; n, m)Q(n, m; n , m!) + (13) 

n-l m+(n— n") 

+ Y, X! K(n,m;n" ,m")Q(n" ,m";n' ,m') = 0, 

n"=0 m"=m—(n—n") 

where 

Q(n, m; ra', m') = ^ / £ ms (A, n)(dp ss/ {X) - d °P ss/ (A)) £v m > (A, n')- (14) 
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Equation (|T3| ) is a two-dimensional analog of the Gelfand-Levitan integral equations for 
finite-difference equation (p. However, this system is not sufficient to determine the 
coefficients K. A supplementary system of equations is obtained upon substituting ( ]T0| ) 
into the complete relation (^) at n' = n and m! = m 

K~ 2 (n, m; n, m) = Q(n, m; n, m) + 
n-l m+(n— n") 

+ E K~ 1 (n,m;n,m)K(n,m;n\m'')Q(n'\m';n,m). (15) 

n"=0 m"=m—(n—n") 

Now let us find connections between the potential coefficients and orthogonalization coef- 
ficients K(n, m; n', m'). To this end, we insert expression (|T0|) for the polynomial functions 

o 

(p ms (X,n) in terms of y? ms (A,n) into the Schrodinger difference equation (|TJ) 



n— 1 m+(n— 1— n') 

"ram 



«W E E K(n-l,m;n',m') <p m > s (\,n') 

n'=0 m'=m— (n— 1— n') 

n+1 m+(n+l— n') 

+a„+i m 51 X! + X > m 5 n '> m ') ^"»'« n ') 

n'=0 m'=m~(n+l—n') 
m+1 n m"+(n-n') 

+ E V mm"(n)J2 E K(n,m";n',m') %> m , a (A,n') 

m"=m~l n'=0 m'=m" — (n—n') 

n m+(n—n') 

= A E E X(n,m;<m')£w s (A,<). (16) 

n'=0 m'=m—(n—n') 

o 

We transform the r.h.s. of the above relation by means of the substitution of A ^P m 's (A, n) 

o o ° 

from the finite difference equation (g) for f ms (A, n) with the known potentials a nm , V mm ' 
(n) , {vnl — m — l,m,m + 1) 

n m+(n-n / ) 

A E E K(n,m;n',m') <P m r s (\,n') = 

n'=0 m'=m— in— n') 
n m+(n—n') 

+ E E K(n,m;n',m')( a n , m ,(f mls (A,n'- 1) 

n'=0 m' =m— (n—n') 

m'+l 

W s (A, n' + 1) + y^'m" (n f ) V m »s (A, n')). (17) 

m"=m'— 1 

Further, we take advantage of the orthogonality relation (|8J) for the matrix functions 

o o 

$ (A,n) orthogonal with the weight matrix p (A). Multiplying expression (|X6|) with its 

o o 

transformed r.h.s. ( |TT| ) by VVm (A, "ft + 1), integrating over d P ss > (A), and summing up 
over the indices s and s', we arrive at the relationship between the potentials a nm , a nm 
and the coefficients K(n,m;n'm f ) 

o K(n,m;n,m) f ^ 

O'n+lm ^n+lr 



K(n + 1, m; n + 1, m) 
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The relations for the coefficients c nm and b nm are established in a similar manner. To 

o 

determine b nm+ \ = V mm+ i[n), eq.fllq) is multiplied by f s 'm+i (A, n) with (|T7|) taken into 
account and integrated with the weight p ss / (A) by using the orthogonality (§). As a 
result, we have 

o K(n,m;n,m) o K (n, m; n — 1, m + 1) 

Onm+l 6nm+l "777 I i I TT a nm+l "777 I i I 7T 

iv (n, m + 1; n, m + 1) K [n, m + 1; n, m + 1) 

if(n+ l,m;n,m+ 1) , . 

_a n+l m ^77 — j ' TT- 

K (n, m + 1; n, m + 1 j 

o 

The relation for c nm = K nm (n) is derived analogously, only (|lq ) is multiplied by f s ' m (A, n) 
o o K(n,m;n — l,m) K(n+ l,m;n, m) 

Cnm C nm -|- &nm 777 \ "n+lm 777 \ ■ (~^/ 

K{n,m;n,m) K{n,m;n,m) 



Substituting ( |Tq) into (|20f) we arrive at 

o o K{n, m; n — 1, m) o K(n + l,m;n,m) 

Cnm C nm -\- 0i nm — — - CL n _\_i m ■ - ■ - r. 

K{n } m;n,m) K(n + 1, m; n + 1, m) 

o ° 

At a nm =a nm = 1, b nm =b n m= 1, the derived generalized expressions turn into more 
simple ones presented in J7|. The two-dimensional finite-difference inverse problem under 
consideration is also a generalization of that 113 1 with the potential coefficients a nm ^ 1, 



bnm 7^ 1, connected nevertheless in a special way. 

In principle, it is easy to formulate the problem of restoring the matrix V mm > (n) in (|2|) 
with all nonzeroth elements, like in a multichannel problem 0. The latter would corre- 
spond to the potential being nonlocal with respect to one of the coordinate variables (in 
our case "m"). If the consideration were made in the polar coordinate system, nonlocality 
with respect to discrete angles would occur. In the case of continuous coordinates, the 
inverse problem for the potential, nonlocal relative to angles, was considered by Kay and 
Moses ia. 



III. BARGMANN DARBOUX TRANSFORMATIONS FOR THE 
TWO-DIMENSIONAL DISCRETE SCHRODINGER EQUATION 

In this section, we describe the algebraic procedure by taking into considaration simple 
kernels Q in the form of a sum of several terms with a factorized coordinate dependence 

P o 

Q(n,m;n',m') = ^ ip {X^,n,m) ip (X^n'.m!) 

P 

= Y,H^rns (A M ,n)7 s (A M )^7 s ,(A M ) f s , m , (A M , n). (22) 

jU=l S gl 

o o 

Here the functions ip (A M , n, m) =ip m (A M , n) are combined as elements of the vector 

o o o o 

^> (A M ,n) = (... ip x (n),ip 2 ( n );---;^m (Xfj,, n ), ---V obtained as a product of the matrix 
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solutions $ (A M , n) taken at eigenvalues A = A M of the reconstructed H and the vector 

r(A M ) 

i[) (\^,n,m) = 2^ Vms (A M ,n)7 s (A At ). 

s 

The elements 7 S form the normalization matrix C(A M ) = r(A M )r'(A At ) with elements 
Css'(A^) = 7s(A Ai )7s'(A At ) corresponding to the bound state ip^in^m) = ^(A^, n, m), = 
l,2,..p. 

Like Q, the ortogonalization kernel K(n, m; n'm') is also presented as a sum of several 
factorized terms. Really, substituting ( [22] ) into the Gelfand Levitan equation (|TB|), we 
obtain 



p n m+(n-n") Q 

" E( E E #(n,m;n",0 ^ (A„, n", m")} V (A M ,n'X)- 

/Lt=l n"=0 m"=m—(n— n") 



(23) 



Noting that the expression in braces is the solution ip^n, m) (pi]) at A = A M of eq.([L]) with 
the desired potentials a nm , b nm and c nm , we immediately receive 



K(n, m;n ,m 



■ u 



(24) 



It is evident now that the new wave functions (p ms (X,n), determined by fllO|) with the 



kernel K taken in the form fl24J), are related to the old ones <j> ms (A, n) by 

p n m+(n-n') Q 

y2 ms (A,n) = -^Vv( n > m ) E E V* (A M , n', rri) V m ' s (A, n'). 

M n'=0 m'=m— (n— n') 



(25) 



In view of ( 24]) for K(n,m;n' ,m') in eqs. (18|), (|i"9|) and (pO), one can immediately 
write expressions for discrete potentials in the closed form 



EjUi^yfom) ^ (A M ,w,m) 

o 

Eu=i + 1, m) ^ (A M , n + 1, m) 



(26) 



bnm+l bnm+1 



+ a 



n,m+l 



— a. 



n+lm ' 



E^=iVy(^V^) ^ (A M ,ra,m) 

o 

E{Ui Vv( n > + 1) V (-V> n, m + 1) 

o 

EjUi m) V> (A M , n - 1, m + 1) 

o 

EjLi m + !) ^ (V n, m + 1) 

o 

EjUi + 1, m) ip (A M , n, m + 1) 

! O 

Eu=i m + 1) V (Ap, n, m + 1) 



(27) 



and 
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o E£=i j^ijh m) (A M , n-l,m) 
E£=iVv( n > m ) i) (A M ,n,m) 
E^=i Qn+im jy«(n + 1, m) j) (A M , w, m) 

o 

££=i Vv( n + i,m)ijj (Xf,, n+l,m) 
The solutions ?/v( n > m ) have to be found from the Gelfand-Levitan equations (0) and 



(p"5l) taking account of (^) and (|24|). Thus, the operator K(n } m;n' \m!) defined by (|24|) 

transforms the solutions y? sm (^) of eq. @ into the solutions tp sm ( n ) of eq.(|l|), determined 
by ( |T0D or ([2"5"D, with the potentials a nm ,b nm and c nm defined by (^), (0) and 



It is not difficult to see from the definitions fll4]) and ( |12"D that the kernels Q and X 
like (|22l ) and (|24l ) can be obtained provided that the spectral weight functions p(A) and 

o 

P (A) for both the sets of potentials coincide except, for instance, p eigenvalues at A = X u . 



This permits one to construct potentials with p new bound states by using (P6|), (|28 
and ( pT| ) or generate the family of spectral-equivalent potentials whose spectra coincide 

o o 

A/i —Xfi and it is only the normalization factors that are different. In the latter 

case Q is taken in the form 

v 

Q(n, m;n , m ) = 2^2^ ^ (A M , n)(C ss /(A^)- C ss / (A M )) W m > (A M , n ) (29) 

fj,=l ss' 

o 

and the above procedure can be used to construct spectral-equivalent operators H and 
H. In spite of a complicated form of the expressions for the potentials (|26|) , ($7\) and 
(p8|), they are simplified for a large set of particular cases. For example, if we deal with 
the free discrete Schrodinger equation as a reference one, fixed by the choice c nm = 0, 
- 1 

Q"nm bnm = 1- 

Let us now consider another simple case when the Hamiltonians differ only by spectral 
data at one bound state. In this case, the summation over [i in all formulae fl25|) - (|28|) 
vanishes. General solutions ip sm (X,n) from (EH) at arbitrary A can be written as 



n m+(n-n') Q 

ip ms (X,n) = -tp(n,m) Y i> { n \ m ') ^m's (A, n'), (30) 

n'=0 m'=m—(n—n') 

o 

where ip (n, m) is a special solution of the Schrodinger equation @ for the value of the 
spectral parameter A = /i and ip(n,m) is the solution of eq.(JI|) with the same eigenvalue 
A = fi. New potentials are expressed in terms of the known old potentials a nm , c nm , and 

o o 

bnm, functions ip (n,m) and functions ip(n,m) that can be determined from the second 
Gelfand-Levitan equation ([15]) 



o ip{n,m) (n,m) 

O-n+lm — "n+lm j l"-U 

ijj(n + 1, m) ip (n + 1, m) 



i/j(n, m) (o ip (n, m) 
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o jj (ra-l,m + l) o jj (n,m) \ 

~r a nm+l 5 a n+lm ~ I • \y^) 

ip (n, m + 1) ip(n + l,m)/ 

and 

o 

if>(n- l,m) o V> (n, m) 

Cum C nm + O nm Q "n+lm • l"'-'/ 

^ (n, m) ip (n + 1, m) 

It should be noted that relationships between potentials and functions can be obtained 
within the Darboux transformation method or factorised method without using formulae 
of the inverse problem. 

Connection between Darboux transformations and inverse problem ones. It is inter- 
esting to note that transformation (|30|) with one bound state corresponds to Darboux 
transformation for the finite-difference equation ([]]). Indeed, let us search for a solution 
(p ms (X, n) of eq.(|I]) with some initially unknown potentials in the form (0). Next it is 
necessary to find conditions for the potentials a nm , b nm and c nm and special functions 
ip(n,m) at which general solutions <p ms (\,n) specified by (|30| ) will satisfy the discrete 
Schrodinger equation ([!]). Substitute (|30|) into ([!]) 

n+1 n+1— n' 

a n+lm tp(n + 1, m) ^ E ^ ( n '' m ') ^ m ' s n ') 

n'=0 m'=m— (n+1— n') 
n-i m+(n— 1— n') 

+a nm ip(n - 1, m) ^ ^ ^ (n', m') ¥ m ' s (A, ra') 

n'=0 m'=m— (n— 1— n') 
-m+1 n m"+(n-n') Q 

+ E KwWt/X, m") ^ E ^ (n',m') <P m > s (A,n') = 

m"=m— 1 n'=0 m'=m" — (n— n') 

n m+(n-n') Q 

= A^(n,m)^ £ V (n',m') £ m < s (AX). (34) 

n'=0 m' =m—(n—n') 



Transform the r.h.s. of (|34[) substituting A ty m 's (A,n) from (||) 

n m+(n-n') Q 

Xtp(n,m) ^ E i> {n',m') <P m 's (\n') = 

n'=0 m'=m— (n— n') 
n m+(n-n') Q 

A^(n, m) 51 E ^ m> ) [^n'm'^m's (A, n' - 1)+ a„' + i m ^ m ' s (A, n + 1) 

n'=0 m'=m—{n—n') 
o o o o o o 

+ <V m /Ws (A,n')+ bn'm'fm'-ls (A,n)+ fe„ m '+lW+l S (A, n )]. (35) 



Further, to obtain the relations for potentials a„+i m , 6 nm+ i and c nm , multiply eq.(|34]) with 

. o o o 

its transformed r.h.s. (pq) by V^ms (A, n + 1), Vm+is (A, n) and ^ ms (A,n) and take into 

O 

considaration the completeness relation (|8]) for the functions <p (A). The expressions for 
a nm , 6 nm and c nm thus derived coincide with formulae (|3l|) , (0) and (|33"D , correspondingly, 
obtained from the formulae of the inverse problem. 
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IV. CONCLUSION 



The Gelfand-Levitan spectral inverse problem for the discrete two-dimensional 
Schrodinger equation is considered on the basis of the Berezanskii technique of orthogo- 
nalization of polynomial matrices. By using the derived formulae of the inverse problem, 
discrete Bargmann-Darboux transformations in two dimensions are given. Analytic rela- 
tionships are established between the solutions for two different sets of discrete potentials 
and the potentials themselves. 

I thank Dr.V.M.Muzafarov and E.P.Velicheva for helpful discussions. 
A part of this work was done in June of 1998 at the Erwin Schrodinger International 
Institute for Mathematical Physics, Wien, Austria. 



Faddeev L.D. 1974 Sov.Probl.Math. 3, 93-180. 

Newton R.G. 1980 J.Math.Phys. 21 1698; 1981 22, 631; ibid, 2191. 

Newton R. 1989 Inverse Schrodinger Scattering in Three Dimensions (Springer- Verlag. 
Berlin Heidelberg/New York). 

Novikov R.G. and Henkin G.M. 1987 Usp.Mat.Nauk 42 93-151 (in russian). 
Berezanskii Yu.M. 1965 The Expansion into the Set of Eigenfunctions of Selfadjoint Oper- 
ators (Naukova Dumka, Kiev) 798 p. 

Spiridonov V. and Zhedanov A. 1995 Methods Appl. Anal. 2, 369-98 

Zakhariev B.N. and Suzko A. A. 1990 Direct and inverse problems. (Potentials in quantum 

scattering) (Springer- Verlag. Berlin Heidelberg/New York, 223p.) 

Suzko A. A. 1985 Physica Scripta 31 447-449; 1986 34 5-7. 

Suzko A. A., 1997 Int. J.Mod.Phys. A, 12 277. 

Matveev V.B. and Salle M.A. 1979 Lett. Matt. Phys. 3 425. 

Matveev V.B. and Salle M.A. 1991 Darboux transformations and Solitons Springer Series 
in Nonlinear Dynamics (Springer- Verlag, Berlin Heidelberg/New York); 
Spiridonov V. 1997 J.Phys.A, Math.Gen. 30 L15-L20. 
Case K.M. and Kac M., 1973 J.Math.Phys., 14 594-603. 

Kay I. and Moses H.E. 1961 JVuovo Chim. 22 689; 1961 Comm.Pure Appl.Math. 14 435. 
Suzko A.A., 1978 Dissertation (1 degree), JINR, Dubna, 106. 



11 



